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Abstract 

Background: Physiologic signals, such as cardiac interbeat intervals, exhibit complex fluctuations. However, 
capturing important dynamical properties, including nonstationarities may not be feasible from conventional time 
series graphical representations. 

Methods: We introduce a simple-to-implement visualisation method, termed dynamical density delay mapping 
("D3-Map" technique,) that provides an animated representation of a system's dynamics. The method is based on 
a generalization of conventional two-dimensional (2D) Poincare plots, which are scatter plots where each data 
point, x(n), in a time series is plotted against the adjacent one, x(n + 1). First, we divide the original time series, x(n) 
(n= 1,..., A/), into a sequence of segments (windows). Next, for each segment, a three-dimensional (3D) Poincare 
surface plot of x(n), x(n+ 7), h[x(n),x(n + 1)] is generated, in which the third dimension, h, represents the relative 
frequency of occurrence of each {x(n),x(n + V) point. This 3D Poincare surface is then chromatised by mapping the 
relative frequency h values onto a colour scheme. We also generate a colourised 2D contour plot from each time 
series segment using the same colourmap scheme as for the 3D Poincare surface. Finally, the original time series 
graph, the colourised 3D Poincare surface plot, and its projection as a colourised 2D contour map for each segment, 
are animated to create the full "D3-Map." 

Results: We first exemplify the D3-Map method using the cardiac interbeat interval time series from a healthy 
subject during sleeping hours. The animations uncover complex dynamical changes, such as transitions between 
states, and the relative amount of time the system spends in each state. We also illustrate the utility of the method 
in detecting hidden temporal patterns in the heart rate dynamics of a patient with atrial fibrillation. The videos, 
as well as the source code, are made publicly available. 

Conclusions: Animations based on density delay maps provide a new way of visualising dynamical properties of 
complex systems not apparent in time series graphs or standard Poincare plot representations. Trainees in a variety 
of fields may find the animations useful as illustrations of fundamental but challenging concepts, such as 
nonstationarity and multistability. For investigators, the method may facilitate data exploration. 
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Background 

Physiologic and physical systems often generate highly 
complex output signals. An inherent feature of these signals 
is nonstationarity, defined as changes in their statistical 
properties (mean, variance, and higher moments) over time 
[1]. Furthermore, nonstationarities are important because 
they may contain information of both basic and trans- 
lational interest. However, conventional representations, 
such as time series graphs, may not fully capture these 
time-varying properties. This challenge motivates the deve- 
lopment of alternative ways to visualise the dynamics of 
complex time series. 

We present a simple, new visualisation approach, based 
on the concept of delay maps, which highlights nonstatio- 
narities and related features. A classical delay map is the 
Poincare plot a [1], widely used in biomedicine by investi- 
gators probing heart rate variability [2]. Our dynamical 
density delay map method, termed D3-Map, extends this 
concept to generate animated, colourised two and three- 
dimensional representations. 

This paper is intended to serve both as a theoretical intro- 
duction and a practical tutorial so that readers can recreate 
all the figures and animations using our data or their own. 
Therefore, all videos, data and source code are made pub- 
licly available (http://reylab.bidmc.harvard.edu/.D3Map/). 

Although Poincare plots are among the simplest rep- 
resentation of a systems phase space [1], they may pro- 
vide important information about the dynamics. For 
illustration purposes here we use heartbeat time series, 
i.e., a sequence of consecutive ventricular interbeat 
(QRS) intervals (termed RR) intervals obtained from 
an electrocardiographic (ECG) recording [2]. However, 
our method is applicable to any time series of sufficient 
length. 



A Poincare plot (map) of a time series is the scatter plot 
of x(n) versus x(n + 1) or equivalently, x(n-l) versus x(n), 
hence the general term delay map° [1-7]. In Figure 1 
(top), we show the RR interval time series from a healthy 
subject during sleeping hours, along with its histogram 
and the standard ("black and white") Poincare plot. As de- 
scribed further below, this dataset was selected because it 
exhibits a common type of nonstationarity characterized 
by relatively abrupt changes of "state". The "comet-like" 
shape of the Poincare plot of the RR interval time series 
(Figure 1, top) indicates a statistical relationship between 
consecutive data points where a given RR interval is likely 
followed or preceded by one of similar duration. In fact, 
the spread of the data points in the direction perpendicu- 
lar to the diagonal of the Poincare plot is a measure of the 
time series' local variance ("short term variability"). 

To highlight the added information that Poincare plots 
provide, we shuffled (randomised) the order of the RR 
intervals. The histogram for this shuffled time series is 
unchanged from that of the original signal. However, the 
Poincare plot is dramatically altered by the breakdown 
of correlations caused by the randomisation procedure 
(Figure 1, bottom). 

This example illustrates how the standard Poincare 
plot resolves one limitation of the histogram, namely that 
the latter does not provide information about correlations 
among the data points. However, standard ("black and 
white") Poincare plots have a salient limitation of their 
own; namely, they do not capture information about 
the density of the data points. This limitation can be 
overcome by modifying the standard Poincare plot to 
include the relative frequency of pairs of consecutive 
data points by constructing a 2D histogram 0 of RR(n), 
RR(n + 1) (Figure 2). 



Time series 



Histogram 



Poincare plot 
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0.6 0.8 1 1.2 1.4 
RR(n), sec 



Figure 1 Time series graph (left), histogram (center) and Poincare plot (right) of the interbeat interval [RR(n)] time series derived from 
an electrocardiographic recording of a healthy subject during sleep (top) and from a randomized/shuffled sequence of these intervals 
(bottom). Note that the histogram of the shuffled time series is identical to that of the original one, while the Poincare plot is substantially 
altered (see text). 
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However, even this so-called "3D Poincare plot" omits 
information about the time evolution of a systems dy- 
namics. For example, consider a time series from a dy- 
namical system with two states, A and B. The 3D 
Poincare plot constructed from the entire time series 
will be the same whether the system spends the first half 
of the time in state A and the second half in state B, or 
whether the system changes states on a much shorter 
time scale. Thus, the 3D Poincare plot for the entire 
time series does not capture the time sequence of state 
changes. This limitation is particularly important when 
probing the dynamics of physiologic systems, which are 
typically non-stationary. Here, we introduce a visualization 
technique, termed D3-Map, based on a moving window 
analysis of the data. 

Methods 

The D3-Map approach consists of five simple steps: 

I. Divide the initial time series into a sequence of either 
overlapping or non-overlapping segments. Here, we 
analyse overlapping segments selected using a 500 sec 
moving window, shifted 10 sec at a time. 

II. For each data segment, construct a normalised 2D 
histogram of RR(n) vs RR(n + 1) and its contour 
map. Let m be the number of counts in the highest 
bin of the histogram. The normalisation procedure 
consists of dividing the number of counts in each 
bin of the histogram by m. Here, we also opted for 
smoothing the histogram, using the dscatter d Matlab 



function with the default parameter, which 
implements the method described in [8]. 
IILColourise both the normalized 2D histogram and its 
contour map. In this implementation, we selected 
the jet colourmap, e where blue represents the lowest 
and dark red-brown the highest density values, 
respectively. The colourised, normalised and 
smoothed 2D histogram of RR(n), RR(n + 1) is called 
the 3D surface Poincare plot. 

IV. In a single figure, present the 3D surface Poincare 
plot, its colourised contour map f and the graph of 
the original data segment used to generate them. 

V. Sequentially present the figures with the analysis of 
each data segment as an animation [5]. Here, we 
used free video editing software VirtualDub (www. 
virtualdub.org) to create AVI movies and animated 
GIF images. g 

The cardiac interbeat interval time series used in the ex- 
amples were from open-access, deidentified data freely 
available at the NIH-sponsored PhysioNet website (www. 
physionet.org) [9]. The links to the specific PhysioNet data 
sets used are indicated in parenthesis in the text below. 

Results and discussion 

The D3-Map visualisation method was first applied to the 
cardiac interbeat interval time series from the healthy sub- 
ject whose data (Additional file 1) are shown in Figures 1 
and 2. All figures presented here for this subject were gen- 
erated from the same data set (nsrOOl at http://physionet. 
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Figure 3 Colourised contour map for the time series presented 
in the top panel of Figure 1. The relative density of each pair of 
points (RR(n), RR(n + 1)) is represented using the jet colourmap. Two 
states, S1 and S2 (local density maxima), separated by a dashed line 
drawn through the region of minimum density, are discernable. 



org/physiobank/database/nsr2db/), which we chose be- 
cause it exhibits a type of nonstationarity characterised by 
relatively abrupt state transitions, which are commonly 
observed in the output of "free-running" physiologic 
signals. 

In this specific case, (Figure 3, colourised 2D contour 
plot), there are two states (termed SI and S2), centered 
around RR intervals of 0.65 sec and 1 sec (corresponding 
to periods of faster and slower heart rates of about 92 
beats/min and 60 beats/min, respectively). These two 
states likely correspond to different sleep stages [10]: the 
faster rates to periods of rapid eye movement (REM) 
sleep (or transient awakenings) and the slower rates to 
deeper sleep. A movie of the 2D contour plots (Anima- 
tion 1, Sleep2D.avi) providing information about the na- 
ture (rapid or gradual) of the transitions between the 
two states is included in the supplemental material. 
From a physiologic point of view, the "flipping" between 
these two states represents a highly nonstationary behav- 
ior, reminiscent of bistability [11]. 

Figure 4 shows "snapshots" of the D3-Map animation 
included in the online supplementary material (Animation 
2, Sleep3D.avi) for four different data segments. These 
snapshots demonstrate system transitions between differ- 
ent states over time. The D3-Map animations show that 
the states, themselves, are not static. Instead they appear 
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Figure 4 Selected snapshots from the D3-Map animation derived from the analysis of the interbeat interval time series shown in the 
top panel of Figure 1. Each snapshot comprises: a 500 sec interbeat interval time series segment (bottom), its colourised 2D contour map 
(middle) and colourised 3D surface Poincare map (top). (See text for details). 
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to consist of multiple sub-states, manifesting as local un- 
dulations in the 3D surface, which continuously emerge 
and disappear. This finding is not apparent from visual in- 
spection of the raw interbeat interval time series. 

To further explore the utility of this visualisation 
method, we also show the analysis of cardiac interbeat 
interval time series from two patients with chronic atrial 
fibrillation, a common cardiac arrhythmia that may be 
associated with stroke or heart failure. The RR dynamics 
during atrial fibrillation are presumed to be random. In 
fact, clinicians use the expression "irregularly irregular" 
to describe the unpredictable timings of the heartbeat 
during this arrhythmia. This assumption is supported by 
the circular shape of the contour map in Figure 5 from 
the first of these subjects (dataset alnn from http:// 
www.physionet.org/challenge/chaos/), which is typical of 
time series of uncorrelated values. The D3-Map movie 
(snapshot shown in Figure 6) shows that in this case, the 
dynamics is relatively stationary. 



In contrast, for the second patient with chronic atrial 
fibrillation (dataset a2nn from http://www.physionet.org/ 
challenge/chaos/), the contour map (Figure 7) shows four 
loci of high density consistent with a more structured pat- 
tern of ventricular response, which are not readily dis- 
cerned from the time series. The structure centered around 
(0.4, 0.4) sec in the contour map indicates that the RR in- 
tervals fluctuate around 0.4 sec for at least three consecu- 
tive beats. Similarly, the structure centered around (0.8, 0.8) 
sec indicates that the RR intervals fluctuate around 0.8 sec 
for at least three consecutive beats. The changes from 
0.4 to 0.8 and from 0.8 to 0.4, respectively, generate the 
structures centered around (0.4, 0.8) and (0.8, 0.4) sec. 
The D3-Map (snapshot shown in Figure 8) animation con- 
firms the dynamical nature of these fluctuation patterns. 

This interesting and anomalous pattern of heartbeat 
timings during AF illustrated in Figures 7 and 8 has 
been noted by other investigators [12]. However, its elec- 
trophysiologic mechanism (possibly involving dual AV 
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Figure 5 Colourised contour map (top) of the RR interval time series (bottom) from a patient with typical pattern of atrial fibrillation. 

The circular shape of the contour map is consistent with uncorrelated dynamics. 
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Figure 6 Snapshot from the D3-Map animation derived from 
the analysis of a 1200 sec RR interval time series segment of 
the patient with a typical atrial fibrillation. Top: colourised 3D 
surface Poincare map; middle: colourised 2D contour map; bottom: 
interbeat interval time series. 



nodal pathways or His bundle alternans) and its bedside 
(diagnostic and prognostic) implications have not been 
fully addressed. (Clinical details are not available in the 
two cases presented here). The D3-Map visualization 
tool may help foster basic and clinical work on different 
dynamical subsets of AF, which in current practice are 
all grouped together. 

Conclusions 

We present a new animation method for visualising 
the temporal evolution of complex signal dynamics. For 
trainees, such simple-to-construct plots (termed "D3- 
Maps") may be used to illustrate technical concepts, such 
as nonstationarity and multistability, which are difficult to 
grasp, but essential to understanding the dynamics of 
physiologic control systems. In addition, the animations 
may reveal unexpected patterns in data structure, which 
make the D3-Maps useful in exploratory research, facili- 
tating hypothesis generation and development and testing 
of mathematical/physiologic models. 

Animations 

All animations are available at http://reylab.bidmc.har- 
vard.edu/.D3Map/. 

Animation 1 (Sleep2D.avi) Movie of the 2D contour 
Poincare plots. Each frame shows both the colourised 2D 
contour plot and the 500 sec interbeat interval time series 
segment used to generate it. Movie frames were generated 
using "RR.dat" dataset (Additional file 1) and "Sleep2D.m" 
Matlab script (Additional file 2). 

Animation 2 (Sleep3D.avi) D3-Map animation. Each 
frame shows a 500 sec interbeat interval time series seg- 
ment (bottom) and the corresponding colourised 3D 
surface Poincare (top) and colourised 2D contour 




Figure 7 Colourised contour map (top) of the RR interval time 
series (bottom) from a patient with an atypical pattern of 
atrial fibrillation. 



(middle) maps. Movie frames were generated using "RR. 
dat" dataset (Additional file 1) and "Sleep3D.m" Matlab 
script (Additional file 3). Note: The duration of the 
Sleep2D.avi and Sleep3D.avi movies is 3 min 20 sec, repre- 
senting approximately 6 hours of heart rate dynamics ob- 
tained during sleep. A 500 sec moving window shifted 10 
sec at a time was used to select the data segments dis- 
played in each frame. 

Animation 3 (Typical_AF_2D.avi) Movie of the 2D con- 
tour Poincare plots for a patient with a typical pattern of 
atrial fibrillation. Each frame shows both the colourised 
2D contour plot and the 1200 sec interbeat interval time 
series segment used to generate it. Movie frames were 
generated using "Typical_AF.dat" dataset (Additional 
file 4) and "Typical_AF_2D.m" Matlab script (Additional 
file 5). 

Animation 4 (Typical_AF_3D.avi) D3-Map animation 
derived from the analysis of the RR interval time series of 
a patient with a typical pattern of atrial fibrillation. Each 
frame shows a 1200 sec interbeat interval time series 
segment (bottom) and the corresponding colourised 3D 
surface Poincare (top) and colourised 2D contour (middle) 
maps. Movie frames were generated using "Typical_AF. 
dat" dataset (Additional file 4) and "Typical_AF_3D.m" 
Matlab script (Additional file 6). 
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Figure 8 Snapshot from the D3-Map animation derived from 

the analysis of a 1200 sec RR interval time series segment of 

the patient with an atypical pattern of atrial fibrillation. Top: 

colourised 3D surface Poincare map; middle: colourised 2D contour 

map; bottom: interbeat interval time series. 
^ J 



Animation 5 (Atypical_AF_3D.avi) Movie of the 2D 
contour Poincare plots for a patient with an atypical pat- 
tern of atrial fibrillation. Each frame shows both the 
colourised 2D contour plot and the 1200 sec interbeat 
interval time series segment used to generate it. Movie 
frames were generated using n Atypical_AF.dat" dataset 
(Additional file 7) and n Atypical_AF_2D.m" Matlab script 
(Additional file 8). 

Animation 6 (Atypical_AF_3D.avi) D3-Map animation 
derived from the analysis of the RR interval time series 
of a patient with an atypical pattern of atrial fibrillation. 
Each frame shows a 1200 sec interbeat interval time series 
segment (bottom) and the corresponding colourised 3D 
surface Poincare (top) and colourised 2D contour (middle) 
maps. Movie frames were generated using n Atypical_AF. 
dat" dataset (Additional file 7) and n Atypical_AF_3D.m" 
Matlab script (Additional file 9). 

Note: The duration of the Typical_AF_2D.avi and Typica- 
l_AF_3D.avi movies is 56 sec, representing approximately 
24 h of heart rate dynamics. The duration of the Atypica- 
l_AF_3D.avi, and Atypical_AF_3D.avi movies is 47 sec, 
representing approximately 20 h of heart rate dynamics. A 
1200 sec moving window shifted 300 sec at a time was used 
to select the data segments displayed in each frame of the 
movies. All Matlab scripts listed above use dscatter2 Matlab 
function (Additional file 10) to calculate smoothed Poincare 
map. 

Endnotes 

a The Poincare plot is sometimes referred to as a Lorenz 
plot, a delay map, or a return map. 

b Matlab code for creating Poincare plot from the se- 
quence of RR intervals is listed below: 



data = loadCRR.daf); % load the Additional file 1 with 
the RR time series 
RR = data(:,2); 
RRx = RR(l:end-l); 
RRy = RR(2:end); 

c Due to the normalisation, the height of the 2D histo- 
gram is not exactly equivalent to the probability density, 
but represents a relative occupancy of the phase space 
mapped onto the [0:1] interval. 

d Available at www.mathworks.com/matlabcentral/fileex- 
change/8430-flow-cytometry-data-reader-and-visualisation. 

e A colourmap is defined as a specific chromatic sequence, 
which is linearly mapped onto the [0:1] interval, so any 
number between 0 and 1 corresponds to a particular colour. 

f Since we cannot see the "back" of the 3D surface Poin- 
care plot, we included its 2D contour Poincare plot in 
each figure. 

g We opted to create animations using stand-alone spe- 
cialised software, rather than Matlab see Appendix in 
ref. [13]. Specialised multimedia editing software has 
additional options for post-processing, such as compres- 
sion algorithms, annotation editing (e.g., adding extra 
frames with a movie title, explanatory text or subtitles), 
as well as audio (e.g., narrations or music). 

Additional files 



Additional file 1: (RR.dat). ASCII file with the cardiac interbeat (RR) 
intervals obtained from the electrocardiographic recording of a healthy 
subject during sleep (~6 hour long). This dataset was used to create 
Sleep2D.avi and Sleep3D.avi movies. 

Additional file 2: (Sleep2D.m). Matlab script used to generate the 
frames of the Sleep2D.avi movie. 

Additional file 3: (Sleep3D.m). Matlab script used to generate the 
frames of the Sleep3D.avi movie. 

Additional file 4: (Typical_AF.dat). ASCII file with the cardiac interbeat 
(RR) intervals obtained from the electrocardiographic recording of a 
subject with atrial fibrillation. 

Additional file 5: (Typical_AF_2D.m). Matlab script used to generate 
the frames of the Typical_AF_2D.avi movie. 

Additional file 6: (Typical_AF_3D.m). Matlab script used to generate 
the frames of the Typical_AF_3D.avi movie. 

Additional file 7: (Atypical_AF.dat). ASCII file with the cardiac 
interbeat (RR) intervals obtained from the electrocardiographic recording 
of a subject with atrial fibrillation, whose contour and D3 maps show 
atypical dynamical patterns. 

Additional file 8: (Atypical_AF_2D.m). Matlab script used to generate 
the frames of the Atypical_AF_2D.avi movie. 

Additional file 9: (Atypical_AF_3D.m). Matlab script used to generate 
the frames of the Atypical_AF_3D.avi movie. 

Additional file 10: (dscatter2.m). Matlab function that creates a scatter 
plot coloured by density. This is a modified version of the original 
dscatter function http://www.mathworks.com/matlabcentral/ 
fileexchange/8430-flow-cytometry-data-reader-and-visualization/content/ 
dscatter.m. All Matlab scripts listed above use this function. 



Abbreviations 

D3-Map: Dynamical density delay map; HR: Heart rate. 
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